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ABSTRACT 

We present measurements of the spin-orbit alignment angle, A, for the hot Jupiter systems WASP- 
32, WASP-38, and HAT-P-27/WASP-40, based on data obtained using the HARPS spectrograph. We 
analyse the Rossiter-McLaughlin effect for all three systems, and also carry out Doppler tomography 
for WASP-32 and WASP-38. We find that WASP-32 (T eS = 6140l?g K) is aligned, with an alignment 
angle of A = 10.5° + 6 6 5 4 obtained through tomography, and that WASP-38 (T cff = 6180t^K) is also 

aligned, with tomographic analysis yielding A = 7.5°j~ 6 ][. This latter result provi des an order of 
magnitude improvement in the uncertainty in A compared to the previous analysis of iSimpson et all 
IpOll . We are only able to loosely constrain the angle for HAT-P-27/WASP-40 (T eS = 5190±^ K) to 
A = 24.2°^J|'5 , owing to the poor signal-to-noise of our data. We consider this result a non-detection 
under a slightly updated version of the alignment test of iBrown et al.l (|2012h . We place our results in 
the context of the full sample of spin-orbit alignment measurements, finding that they provide further 
support for previously established trends. 

Subject headings: planets and satellites: dynamical evolution and stability - stars:individual(WASP- 
32, WASP-38, WASP-40) - techniques: radial velocities - techniques: spectroscopic 



1. INTRODUCTION 

The number of planetary systems for which the sky- 
projected, spin-orbit alignment angle has been measured 
is steadily increasing, and is now approaching the point 
at which serious statistical analyses can be made. The 
majority of these angles have been measured through the 
Rossiter-McLaughlin (RM ) effect flHoltl|1893t [Schlesingerl 
119101 119161 : lRossiterlll924t lMcLaughlinlll924D . a well es- 
tablished technique that considers the small anomaly in 
the radial vel ocity curve that is produced by a tran- 
sit event fe.g IQueloz et al.ll2000a[) . However there are 
a growing number of systems for which the misalign- 
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ment angle has been measured using alternative means. 
In some cases this is out of necessity, whilst in oth- 
ers it arises from a desire to expand the repertoire of 
analysis methods that are available, in an effort to re- 
duce the ever-increasing demands on spectroscopic in- 
struments. Examples of the alternatives currently avail- 
able include analysis of the effect of s tar spots on the pho- 
tome tric transit observations (e.g. iSanchis-Oieda et all 
120111), consideration of the effect of gravity darkening 
(e.g. iBarnes et ail 1201 ID . comp arison of the mea sured 
and predicted stellar v sin I fe.g.lSchlaufmarj[20T0l), an d 
Doppler tomography (e.g. lCollier Cameron et al.ll2010ap . 
This last method, whilst not greatly reducing the tele- 
scope time required, is able to break the degeneracy be- 
tween the sky-projected alignment angle and stellar ro- 
tation velocity in systems with low impact parameter. 
It is best suited to analysing hot, rapidly rotating exo- 
planet host stars, although it can be applied to planetary 
systems with a range of host parameters. In fact all of 
the alternative methods are complementary to the RM 
measurement approach, allowing as they do the study of 
systems with vastly different properties, and with which 
that traditional method struggles to cope. 

It is becoming increasingly important to push the 
boundaries of the explored parameter space in this way. 
The spin-orbit alignment is an excellent diagnostic for 
competing theories of planetary system formation and 
exoplanet migration; as the number of systems for which 
it is measured increases, so too does our understanding 
of these processes. 

The generally accepted scenario has hot Jupiters form- 
ing beyond the 'snow line ' and migrating inward s to 
their observed separations ( Sassclo v fc Lecarll2000D . It 
is the process by which this migration occurs that is 
disputed. Loss of angular momentum through interac- 
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tions with a protoplanetary disk (|Goldreich fc Tremainel 
1980) was initially proposed as the dominant mechanism, 
which, if it is assumed that such disks are well-aligned 
with the stellar spin axis, would produce a population of 
hot Jupiters in well-aligned orbits. It is worth pointing 



out, however, that this assumption of aligned protoplane- 
being chg 
peers. Lh 

vestigated (|Watson et al.ll201lD . The discovery of hot 



tary d isks is increasin g ly being challenged (e.glBate et al.l 
[20101: lLai et al.l I20T1 iRoeers. Lin fe Laul 120121 ) and in- 



Jupiters in strongly misaligned orbits, including some 
that are orbiting in a retrograde direction, has also led 
to the development of comp eting theorie s that utilise 
the Kozai-Lidov mechanism (iKozail 119621: iLidovl 119621 : 



iFabrvckv fc Tremainel l2007t iNaoz et al.l 1201 ID planet 
planet scattering (fWeiden schillin g fe Marzarj I1996D . 
tidal fr iction, or some com b inatio n of these pro- 
cesses (IFabrvckv fc Tremainel 120071 iNagasawa et al.1 

' >012h. 



2008; INaoz. Farr fc Rasiol l2012). These mechanisms are 
naively expected to produce spin-orbit alignment distri- 
butions that are closer to isotropic. However the true 
picture has turned out to be more complex, and appears 
to lie somewhere between these two extreme distribu- 
tions; 

|W inn et all (j2010al ) found an apparent link between 
stellar effective temperature and alignment angle; plan- 
ets in misaligned orbits seem to preferentially orbit 'hot' 
stars (T c ff > 6250 K), whilst aligned planetary orbits 
seem to be found mostly around 'cool' stars (T c s < 
6250 K). They suggested that this might be connected 
to the size of the convective envelope, with tidal realign- 
ment of orbits around 'hot' stars being s uppressed ow ing 
to their small convective zone. This led |W inn et all to 
conclude that the A distribution (at the time of their 
publication) could be explained by the combination of 
planet-planet scattering and Kozai-Lidov cycles. 

Another potential pattern in the data, and one that 
was identified early in the development of this sub-field, 
was that planets with high mass tend to b e misaligned 
(but not retrograde) ((Johnson et al.l 120091 ) . Counter- 
examples to the initial trend have been found (for ex- 
ample HAT-P- 7 dWinn et al.1 120091: iNarita et all 12009ft 
and WASP-18 (jHellier et all 120091: iTriaud et all ' " 



but planets with M p > 3Mj up do seem to have a 
different distribution of spin-orbit misalignment angles 
(jHebrard et al.H20111l) . This was tentatively interpreted 
as possible evidence for a combination of Kozai-Lidov 
cycles with tidal circularisation and realignment, but 
small number statistics were cited as a cautionary fac- 
tor. More recent work along similar lines has tended to 
concentrate on the ratio of the planetary mass to the 
stellar mass rather th an the planetary mass in isolation 
(jAlbrecht et alJl2012t ). 

A more recently discovered cor relation i s that of align- 
ment angle with host star age. ITriaud! (j2011l ) noticed 
that, for stars with A/* > 1.2A/©, all systems older 
than 2.5 Gyr are well-aligned. This implies that the dis- 
tribution of A changes with time, which in turn sug- 
gests that some misalignment mechanism must oper- 
ate during the youth of hot Jupiter systems, followed 
by some method of realigning the system as it evolves. 
If age is the primary factor then having tidal interac- 
tions as the governing mechanism for this latter stage 
would fit with the observed age trend, as planets around 
older stars will have had longer to tidally realign. On 



the other hand it may also be that strongly misaligned 
planets are simply being destroyed much more quickly 
than their aligned cousins; indeed, such an effect has 
been theoretically demonstrated for retrograde planets, 
which are predicted to reach disrupti on distances sev- 
eral times faster than pro grade planets ( Bark er fc Ogilviel 
120091: iWinn et alJl2010aT ). This would lead to a d ecrease 
in the number of hot Jupiters with time, yet iTriaudl 
(|2011[) found no such trend. Either tidal realignment 
occurs faster than orbital decay, or some other mecha- 
nism is responsible for the evolution of the distribution 
of angles that we observ e. 

Albr echt et al.l (|2012ft have re-examined all of these 
previously detected trends using an updated, homoge- 
neous database of Rossiter-McLaughlin measurements 
that included their own new measurements and re- 
analyses. They found that all of the existing trends are 
consistent with the idea that tidal interactions are re- 
sponsible for the evolution of the spin-orbit alignment in 
hot Jupiter systems. They also considered the depen- 
dence of A on the scaled orbital distance, finding that it 
too is consistent with a tide-driven evolutionary picture. 
Their estimates of characteristic tidal timescales showed 
that systems which were expected to align rapidly exhibit 
angles consistent with alignment, whilst those for which 
tidal realignment was predicted to be weaker display 
a nearly random distribution of angles. Al brecht et all 
stop short however of claiming any mechanism for the 
production of the initial distribution of A which, from 
evidence collected so far, seems to be required to be 
isotropic. 

Despite all of this, the question of how hot Jupiters 
appear where they are and with the spin-orbit angles 
that they have is far from settled. There is only so much 
that 'typical' transiting hot Jupiters can tell us; it is 
the more unusual systems, lying at the extremes of the 
distributions in mass, effective temperature, and wsinJ, 
and that can only be accessed through methods such as 
Doppler tomography, that will provide the best test of 
the theory underlying the evolution of orbital misalign- 
ment with time. It is also important to use newer anal- 
ysis methods to examine systems in tandem with the 
consideration of the RM effect in order to get to grips 
with their intricacies, strengths, weaknesses, and inher- 
ent error characteristics. In this paper we present new 
measurements of the spin-orbit alignment angle for three 
WASP systems. Two, WASP-32 and HAT-P-27/WASP- 
40 (hereafter WASP-40), have not previously had been 
analysed. The third system, WASP-38, has been ex- 
amined before using the RM effect, but we present new 
spectroscopic data that improves upon the existing pa- 
rameter uncertainties. For all three systems we compare 
the results obtained using the RM effect and Doppler 
tomography. Throughout we will characterise the spin- 
orbit alignm ent angle as A, foll owing the convention es- 
tablished by lOhta et all (|2005l) and widely followed in 
the literature, r ather than the alter native convention of 
/3 = -X used bv ITriaud etHI (pOlOf l. 

2. DATA ANALYSIS 

2.1. Rossiter-McLaughlin measurements 

Our analysis method is base d around an adapted 
version of the code described in [Collier Cam eron et all 



The alignment of WASP-32, WASP-38 and WASP-40 
TABLE 1: Existing system parameters for WASP-32, WASP-38, and WASP-40. 
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Parameter 


Unit 


WASP-32 


WASP-38 


WASP-40 


M, 


Mq 


1.10 ±0.03 


1.203 ±0.036 


0.921 ±0.034 


R* 


Rq 


1.11 ±0.05 


1 'iQI +0.030 
±.ooi_ 025 


0.64 ±0.031 


Teff 


K 


6100 ± 100 


6150 ± 80 


5246 ± 153 


v sin I a 


km s 


K K 4- n A 
0.0 It U.4 


o.o it U.4 


z.4 it U.O 


~, b 




o.o in u.o 


O. 1 m u.o 


i n -t- n ^ 

1 .U It u.o 


M p 


MJup 


3.60 ±0.07 


2.691 ±0.058 


0.617 ±0.088 


Rp 


i?,Jup, eq 


1.18 ±0.07 


1 094+ 029 
i.uy^_ 02g 


l-U3S_o.o50 


P 


days 


2.718659 ± 0.000008 


6.871814 ± 0.000045 


3.0395589 ± 0.0000090 


a 


AU 


0.0394 ± 0.0003 


07522+ - 00074 

U.U/OZZ_ Q 00075 


0.03995 ± 0.00050 


e 




0.018 ± 0.0065 


0314+ 0046 
u.uo±<t_ 0041 


O(adopted) 


i 


deg 


85.3 ±0.5 


88.83lS;g 


85.0lt° : » 


ref 




1 


2, 3 
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References. — I'll IMaxted et al.ll20TOl : C2'l lBarros etal1l20lTl ; f3l ISimpson et alJI20TJ : (4) I Anderson et al.ll20Tll 

a v sin I have been updated through spectroscopic analysis of the 
new HARPS data. 

%ac values were obtained using the lBruntt et al.l l2010l'l calibra- 
tion against T e g. 



(2007). It utilises a Markov Chain Monte Carlo 
(MCMC) algorithm , and has previously been detailed in 
IBrown et al.l p012). We model the complete sets of pho- 
tometric and spectroscopic data to maintain consistency, 
account for parameter correlations, and fully characterise 
the uncertain ties in our results. Since the publication of 
IBrown et al.l we have made some small updates to the 
code to increase its functionality. Rather than using a 
global stellar 'jitter', individual radial velocity (RV) data 
sets are now allocated 'jitter' values indiv idually. Sim- 
ilarly, the line dispersions required for the iHirano et al.l 
(2011) formulation for modelling the RM effect are now 
calculated for each separate set of RV data. We have also 
up dated our stellar r adius prior to use th e calibration 
of ISouthworthl (|2011l ) rather than that of lEnoch et al.1 
(|2010D . as the former uses a greater number of stars and 
focuses on a mass range that is directly relevant to exo- 
planeta ry systems. 

As in IBrown et al.l (|2012f ). we apply four Bayesian pri- 
ors in all possible combinations in an attempt to fully 
characterise the systems under consideration. We apply 
priors on orbital eccentricity, spectroscopic usin/, long- 
term RV trend, an d stellar radius (using the method of 
lEnoch et al.l (l2010D in conjunct ion with the updated coef- 
ficients fromlSouthworth (2011)). To distinguish between 
the combinations of priors we consider the reduced spec- 
troscopic x 2 , which we refer to as x 2 cd . If there is no 
combination of priors with a significantly lower value of 
Xrcd we choose the model with the fewest free parameters. 
The application of the stellar radius prior we consider on 
the basis of the statistical parameter S (the stellar radius 
penalty) (jCollier Cameron et aTll2007D . 



S = -2]nP(M*,R*) 



(i?„ — Rq)' 



(1) 



where M* and i?* are the stellar mass and radius as 
calculated by the MCMC algorithm, Rq is the stellar 
radius derived from the (J-H) colour, and an is the la 
error in Rq. S measures the discrepancy between the 



two stellar radius values, and if we find a large increase 
in S when the stellar radius prior is removed, we choose 
a solution in which it is applied as our preferred one. 

RV measurements for our new HARPS data were cal- 
culated through a Gaussian fit to the cross-correlation 
functions (CCFs), using a window of three times the 
FWHM. HARPS spectra cover the wavelength range 
378nm < A < 691nm. 

2.2. Doppler tomography 

Our Doppler tomography method also uses the com- 
plete set of photometric and spectroscopic data for 
an exoplanet system, and is again based around a 
modified version of t he M CMC code discussed by 
iCollier Cameron et al.1 (|2007h . In this case however, the 
alignment of the system is analysed through a compari- 
son of the in-transit CCFs (covering the same 378nm < 
A < 691nm range as those used for the RV calculations) 
with a model of the average stellar line profile. This latter 
model is created by the convolution of a limb-darkened 
stellar rotation profile, a Gaussian representing the lo- 
cal intrinsic line profile, and a term corresponding to the 
effect on the line profile of the 'shadow' created as the 
planet transits its host star. This 'bump' is time- variable, 
and moves through the stellar line profile as the planet 
moves from transit ingress to transit egress. The pre- 
cise trajectory of the bump is dictated by the impact 
parameter, b, and spin-orbit alignment angle, A, which 
together determine the precise value for the stellar radial 
velocity beneath the planetary 'shadow' at any moment 
during the transit. This leads to a more accurate model 
of the spectroscopic transit signature than provided by 
RM analysis. 

As noted previously this provides a powerful method of 
analysing spin-orbit alignment that is able to explore pa- 
rameter space unaccessible to the Rossiter-McLaughlin 
method (such as rapidly rotating host stars) whilst 
breaking degeneracies inherent in the other method. It 
has already been used to provide new constraints on 
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the spin-orbit a ngles of the WASP-3 (iMiller et alll201Cj) 
and HD189733 ([Collier Cameron et al. 2010a) systems , 
to analyse the WASP-33 dCollier Cameron et al.ll2010bl) 
and CoRoT-11 b iGandolfi et al.1 (|2012l ) systems? and to 
examine the alignment in five further systems (WASP- 
16, -17, -18, -23 and -31) (Miller et al., accepted). 

3. WASP-32 

WASP-32b is a dense hot Jupiter in a 2.72 day orbit 
around a Sun-like (spectral type G, T G g = 6140t^oo^)' 
lithium depleted star, and is one of only a small number 
of hot Jupiters with a mass gre ater than 3 Jupiter m asses. 
Its discovery was presented bv lMaxted et al.l (120101), who 
used photometry from WASP-S (Poll acco et al.l I2006T) 
and Faulkes Telescope North (FTN), in concert with 
spectro scopic observations f rom the CORAL IE spectro- 
graph (jQueloz et al.l l2000at iPepe et al.l 12002ft . to deter- 
mine the existence of the transiting planet. 

We used the HARPS spectrograph to observe the tran- 
sit of WASP-32 b on the night of 2011 September 26. 
Thirty observations were acquired over the duration of 
the night, and additional data were collected on the 
nights of 2011 September 24, 25 and 27 (see journal of ob- 
servations, Tabled]). We obtained simultaneous photom- 
etry of the same transit using EulerCam, mo unted on the 
1.2m Leonard Euler telescope at La Silla (jLendl et al.l 
[20121) . and using the TRAPPIST telescope at La Silla 
( Jehin et al.l 1201 lT ) . We also obtained photometry of a 
further transit using TRAPPIST, on 2011 November 24. 

We carried out a spectroscopic analysis of the new 
HARPS spectra to determine an updated estimate of 
vsml for the host star. We assumed a macroturbu- 
lence of ^ mac = 3.5 ± 0.3kms _1 using the calibration of 
iBruntt et al.l ([2010D . They describe an analytical poly- 
nomial correlation between T c g and v ma c (their equation 
(9)) by convolving synthetic line profiles with different 
wsin/ and i> ma c values, and fitting to high signal-to-noise 
spectra from several instruments, including HARPS. We 
obtained vsml = 5.5±0.4kms _1 , in agreement w i th the 
value of 4.8 ± 0.8 km s" 1 found by IMaxted et~aT1 pOlOh 
from their CORALIE spectra. Our new value was ap- 
plied as the spectroscopic prior. 

3.1. Rossiter- McLaughlin analysis 

We initially applied a stellar 'jitter' of 1.0 m s _1 to both 
the existing CORALIE RV data and our new HARPS 
data; this value is below the level of precision of the spec- 
trographs used for this work, and was added in quadra- 
ture to any data points falling within the transit. The 
values of x^ed that were returned by our algorithm with 
this level of 'jitter' applied all fall with la of 1.0, indi- 
cating that the solution is well constrained. There was 
therefore no need to increase the level of stellar activity 
accounted for by our modelling. 

Adding a long-term, linear RV trend produced no dis- 
cernible effect on the quality of fit that we obtained, or 
on the value of Xred ■ Relaxing the prior on the stellar ra- 
dius led to only marginal changes in the values of S, M*, 
p„, and b. It also produced no change in the value 
of Xredi we therefore conclude that any such trend is in- 
significant, and choose as our definitive solution a model 
which does not apply the prior. Similarly, we elect not 
to apply the prior on v sin / in our final solution. Whilst 



applying the prior produced an increase of 1.0 km s -1 in 
the value of v sin / returned by the MCMC algorithm, it 
had no impact on the value of X? e d- 

The pr ior on orbital eccentr icity required more careful 
analysis. Maxted et al. (2010) reported a marginal 2.8cr 
detection of eccentricity in the planet's orbit, and sug- 
gested that it could be confirmed through observations of 
the secondary eclipse. To our knowledge no such observa- 
tions have been carried out, so we approach the question 
with all options available. Our models with floating ec- 
centricity all find e < 0.014, s lightly less tha n the value 
of e = 0.018 ± 0.0065 found bv lMaxted et all (poTol). and 
none show any improvement in Xred compared to the 
equivalent models with fixed, circular orbits. We tested 
the significance of the eccentrici ty values recovered by ou r 
algorithm using equation 27 of ILucv fc Sweeney! (| 1 9 7 If) . 
which adopts a null hypothesis of a circular orbit and 
considers an orbit to be eccentric if this is rejected at the 
5 percent significance level. This F-test indicated that 
none of the eccentricities are significant, and thus that a 
circular orbit is favoured. 

Our adopted model thus uses the combination of a cir- 
cular orbit and no long-term RV trend, with neither the 
vsml or stellar radius priors applied. This model pro- 
vides values of A = 8.6°j^' 5 , vsml = 3.9 ± 0.5 km s -1 
(slightly slower than the value from spectroscopic anal- 
ysis), b = 0.66 ± 0.02 and i = 85.1° ± 0.2. The result- 
ing RV curve is displayed in Fig. Q] alongside a close-up 
of the transit region, showing the Rossiter-McLaughlin 
anomaly. The amplitude of the anomaly is low ow- 
ing to the moderate rotation speed of the host star, 
but the signal-to-noise is high and the anomaly is well 
constrained. We found that the semi-amplitudes re- 
turned for all three of th e RV d ata sets (the CORALIE 
data from IMaxted et all (|2010l ). our new HARPS out- 
of-transit data and our HARPS in-transit data) were in 
good agreement, and consistent with the results from the 
discovery paper. Our barycentric velocities on the other 
hand whilst consistent w ith each other , are sl ightly less 
than the value found by IMaxted et al.l (|2010l ) , even for 
their CORALIE spectroscopy. 

3.2. Doppler tomography 

The set of priors identified as comprising the best- 
fitting model for our RM analysis were applied to our 
Doppler tomography method, allowing us to assess a 
single model only. Fig. [2] displays the residual maps 
from our analysis, which returned values of v sin / = 
3.9±g;£ km s" 1 and A = 10.5° + 5 6 9 4 . 

There is little to choose between the results returned 
by our two analysis methods. In this case, since the con- 
straints on the spin-orbit angle were well-defined by our 
original RM analysis, the tomography method has been 
unable to provide much improvement. However it does 
confirm the results from the traditional, RV measure- 
ment based RM analysis, namely that the system is well 
aligned, and in a prograde orbit. This is easily seen in 
Fig. [2] Fig. [2a] shows the time-series map of the residu- 
als after the subtraction of the stellar line profile only; 
the effect of the planet therefore shows up as a bright 
'streak' across the figure, centred on phase and the 
barycentric radial velocity of the host star, and travel- 
ling between the ±usin/ values. The trajectory of the 
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Orbital phase 



Orbital phase 



Fig. 1. — : Results from our adopted model for WASP-32: e = 0; no long-term radial velocity trend; no prior on 
the spectroscopic usin/, and no stellar radius prior. The best-fitting model is plotted as a solid black line. Top left: 
Complete radial velocity reflex motion curve. Data from CORALIE are denoted by triangles. Data from HARPS 
are denoted by squares. Error bars are marked in grey; some are smaller than the size of the data points that they 
accompany. Bottom left: Residuals from the RV fit, exhibiting no correlation with phase. Top right: Close up of the 
transit region from the radial velocity curve showing the RM effect, along with the residuals. Middle right: Close up of 
the transit region, with the orbital contribution removed. Bottom right: Residuals for the radial velocity data within 
the RM window. Colour versions of these figures are available in the online edition of the journal. 



planet signature unambiguously identifies the planetary 
orbit as prograde, moving as it does from bottom-left 
(— vsml at the orbital phase corresponding to ingress) 
to top-right (+wsin/ at the orbital phase corresponding 
to egress). Fig.[2b]in turn displays the final residual map, 
after the removal of the planet signature. The lack of any 
notable, consistent deviation from the mean value of the 
map indicates a lack of significant stellar activity in the 
host star, as any such activity would produce signatures 
similar to that of the planet (e.g. non-ra dial pulsation in 
WASP-33. [Collier Cameron et al]l2010bh . 

4. WASP-38 

The WASP-38 system consists of a massive (2.7Mj up ) 
hot Jupiter in a long (6.87 d), eccentric orbit around a 
bright (V=9.4), rapidly rotating star of spectral type 
F8 and T cB = 6180tggK. Fur ther informat ion regarding 
its discovery can be found in iBarros et al.l fj201 lh . Pho- 
tometry from the WASP-N array, the RISE instrument 
mounted on the 2 m Li verpool Telescope ()Steele et al.l 
I2008t iGibson et al.ll200l ) and an 18 cm Takahashi astro- 
graph at La Palma were combined with spectroscopic 



measurements taken using the CORALIE and SOPHIE 
instruments to confirm the presence of the planet. 

The Rossiter-McLaughlin effect of WASP-38 b has 
been analysed previously by ISimpson et al.l (|2011D , who 
obtained spectroscopic observations of a transit event us- 
ing the FIES spectrograph mounted on the Nordic Opti- 
cal Telescope (NOT) at la Palma. Despite the low pre- 
cision of their measurements, they were able to place 
useful constraints on the misalignment angle using the 
shape of the RV anomaly during transit, ruling out 
high angles and confining the system to prograde orbits. 
They reported a final value for the misalignment angle 
of A = 15°^!!, were not able to provide a firm con- 
clusions as to the alignment, or otherwise, of the system. 

We obtained new spectroscopic observations using 
HARPS of the transit event on the night of 2011 June 
15, as well additional observations made throughout 
2011 to provide coverage of the entire radial veloc- 
ity curve (see journal of observations, Tables^ and [7|) . 
We again obtained photometric observations of a tran- 
sit using TRAPPIST, on 2011 April 13, but unfortu- 
nately we were not able to obtain simultaneous pho- 
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Barycentric RV [km s 1 ] Barycentric RV [km s 1 ] 

Fig. 2. — : Left: Residual map of WASP-32 time series CCFs with the model stellar spectrum subtracted. The 
signature of the planet moves from bottom-left to top-right, supporting the aligned, prograde orbit conclusion from 
our RM analysis. Right: The best-fitting model for the time-variable planet feature has been subtracted, leaving the 
overall residual map. The lack of any features in this figure indicate a lack of large-scale stellar activity. 
The horizontal dotted line marks the mid-transit phase. The vertical dotted line denotes the stellar radial velocity, 
whilst the vertical dashed lines indicate ±wsin7 from this, effectively marking the position of the stellar limbs. The 
crosses mark the four contact points for the planetary transit. 



tometry of our spectroscopically observed event. As 
with WASP-32, we analysed our new HARPS spec- 
tra to obtain a value for wsini of 8.3 ± 0.4kms _1 . 
A macroturbulence of v mac — 3.7 ± O Skms" 1 was 
assumed, again using the calibration of IBruntt et all 
Our v sin I is in excellent agreemen t with the va l ues o f 
v sin I = 8.6 ± 0.4 km s^ 1 quoted by IBarros et all (|2011f) 
and v sin I = 8.58 ± 0.39 km s -1 found by Simpson ct al. 
(|2011|) . Our adopte d v mac is significantl y lower than the 
4.9 ± 0.4 km s " 1 that IBarros et all (|2011D used to fit their 
spectro scopy. IBarros et al.l use cT the calibrat i on of IGravl 
(2008), whereas we used that of IBruntt et"all (|2010D . Re- 
analysing our new spectra using the IGravl calibration re- 
turns a slightly lower value of v sin I = 7.9 ± 0.4 km s -1 , 
in agreement with the IBarros et al.l results. In spite of 
this, we feel that the IBruntt et al.l calibration gives a 
better fit to our data, and it is therefore that result that 
we use for our spectroscopic prior. 

4.1. Rossiter- McLaughlin analysis 

Our initial stellar jitter estimate of 1ms" 1 led to 
poorly constrained results, with the lowest Xied va hie re- 
turned by any of the models being 1.7 . We t h erefor e 
recalculated the stellar jitter following Wright! (2005), 
obtaining three distinct values. We found that in or- 
der to force Xrcd ~ 1 we had to apply the conservative, 
20th percentile estimate of 2.1 ms _1 to our new HARPS 
data, and the 80th percentile estimate of 6.6ms" 1 to the 
pre-existing FIES, CORALIE, and SOPHIE data 11 . 

We found that applying the spectroscopic prior on 
v sin L made little difference to the quality of fit that we 
obtained, or to the values of vs'mL and A that we ob- 
tained when compared to the equivalent case without the 
application of the prior. Similarly, applying a long-term 

11 For an explanation of these estimates, sec Wright (2005) 



RV trend had no effect on the results, and the magni- 
tude of any possible trend was found to be insignificant 
at 7I < 22ms _1 yr _1 . The stellar radius prior however, 
despite producing only a small change in the values of 
Xred> na, d a significant impact on the results that we ob- 
tained. Cases in which the prior was not applied saw av- 
erage increases in the stellar mass and radius of 6 percent 
and 27 percent respectively over their equivalent cases in 
which the prior was applied, as well as an average de- 
crease in the stellar density of 49 percent. Note that these 
changes do not necessarily match perfectly, as under our 
model the stellar density is calculated directly from the 
transit light curves, and independently from the stellar 
mass and radius. Relaxing the prior also produced sig- 
nificant increases in vsml, and substantially raised the 
impact parameter from b = O.I5+030 to b — 0.62+q 1 . 1 ). 
Furthermore, we found that removing the prior increased 
the_value of the stellar radius penalty, S, from S = 14.0 
to S = 105.3. 

Allowing the eccentricity to float led to a clear and 
significant difference in both Xrcd and the total \ 2 f° r 
the combined photometric and spectroscopic model. The 
values returned by our algorithm lie at > 7a from 
e = 0, and were found t o be significant by the F-test 
of lLucv fc Sweenevl (|1971 | ). This co nfirms the eccentric- 
itv detection of IBarros et al.l (|201lD . and the values that 
we find are consistent with the value of e = 0.0314+q 0041 
reported by those authors. 

Our adopted model for this system therefore uses an ec- 
centric orbit, does not include a long-term RV trend, does 
not apply a prior on v sin /, and does utilise a prior on the 
stellar mass. This model returns values of A = 9.2° ti® 5, 
vsml = 7. 7^4 km s -1 (slower than the spectroscopic 

result), b = OMto ll and i = 89.6 + O °|, all of which in- 
dicate a well-aligned system. The radial velocity curves 
are displayed in Fig.[3j The difference in quality between 
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TABLE 2: Comparison of results for WASP-38. 





v sin / 


A 


6 


Source 


(kms- 1 ) 


(deg) 


(•R.) 


Simpson ct al. (2011) 
This work: RM effect 
This work: Tomography 


8.58 ±0.39 

7 7+O.5 
'•'-0.4 
7 c+0.1 
'•°-0.2 


15+ 33 

i0 -43 
q 0+18.1 
a - z -15.5 

7 5+ 4 - 7 
' -°-6.1 


27+°- 10 
u - z ' -0.14 

09+ 13 
12+ 008 



the FIES data presented by iSimpson et al.l (|201lD and 
our new HARPS measurements is immediately appar- 
ent, particularly during the first half of the transit. The 
shape of the anomaly is well defined, and it has the large 
amplitude that is expected given the host star's rapid 
rotation. It also appears to be highly symmetric, lend- 
ing credence to the conclusion that the system is likely 
well-aligned. 

We find that the radial velocity semi-amplitudes 
and barycentric velocities vary somewhat between the 
5 different spectroscopic data sets that we analysed 
(CORALIE data, SOPHIE data, FIES data, HARPS 
data out-of-transit, and HARPS data in-t ransit). In par- 
ticula r, the data obtained using FIES by ISimpson et all 
(|2011h has a much smaller semi-amplitude than any of 
the other data sets; 0.152 ± 0.030 km s -1 compared to 
values betwe en 0.246 ± 0.00 1 and 0.255 ± 0.007kms -1 . 
Interestingly, ISimpson et al.1 found a semi-amplitude of 
0.2538 ± 0.0035 km s -1 in their analysis, but we sus- 
pect that this was overwhelmingly derived from the SO- 
PHIE and CORALIE data, which cover the entire orbital 
phase. The barycentric velocities agree well with the re- 
sults of that previous study however. 

4.2. Doppler tomography 

We again used the set of priors adopted for our RM 
modelling as the basis for our Doppler tomography anal- 
ysis, and Tabled displays the results from th i s ana lysis, 
together with the results from ISimpson et ail (|201lD and 
our own RM analysis. It is immediately apparent that 
we have been able to dramatically reduce the uncertain- 
ties on the projected spin-orbit alignment angle; we will 
return to the question of why this is in Section[7] The 
signature of the planet is clearly defined in Fig. 01 and in 
the final residual image there is no sign of any anomalies 
in the stellar line profiles, indicating that the host star is 
chromospherically quiet. 

5. HAT-P-27/WASP-40 

HAT-P-27 (jBekv et al.l l20lTl) is a fairly typical hot 
Jupiter system, with a 0.6Mj up planet in a 3.04 d or- 
bit around a late-G/early-K type star with T c ff = 
5190j^ 7 Q K and super-Solar metallicity. The system was 
characterised using photometry from HATnet and Ke- 
plerCam on the 1.2m FLWO telescope, and spectroscopy 
from HIRES. It was also independently discovered by the 
WASP survey using the combined WASP-N and WASP- 
S arrays, together wit h spectroscopy from S OPHIE, and 
designated WASP-40 (jAnderson et al.ll2011[ ). 

We obtained new spectroscopic measurements using 
HARPS of the transit on the night of 2011 May 12, and 
carried out additional observations at a range of orbital 
phases throughout 2011 May. New photometric observa- 



tions were also made using TRAPPIST on 2011 May 17, 
covering a full tr ansit. We combine d thes e new data with 
that f rom both IBekv et"aLl (|2011[) and lAnderson et al.l 
(|2011[ ) for our attempt to characterise the RM effect. 

As for WASP-38, we found that our original esti- 
mate of lms -1 for the stellar jitter produced poorly 
constrained (Xr C d ~ models. We calculated possi- 
ble values of 5.9ms -1 (20th percentile), 7.5ms -1 (me- 
di an), and 9.8m s" 1 (80th percentile) using the method 
of iWrightl (|2005l ). and apply the lat ter to the existing 
SOPHIE data. We also note that IBekv et al.l f)2011h 
applied a jitter of 6.3 ms -1 to their HIRES data. To 
confirm that this was reasonable, we analysed the photo- 
metric data in conjunction with only the HIRES spectro- 
scopic data, finding that our initial estimate of 1ms" 1 
produced x^cd = 10.0 ± 1.5, the 20th percentile value 
produced Xrcd = 1-1 ^ 0.5, the median value produced 
Xrcd = 0-7 ^ 0-4, and the 80th percentile value produced 
Xrod = 0.4 ±0.3, whilst applying their estimate produced 
X^d = 1-0 ± 0.5. We therefore follow IBekv et al.1 and 
apply a jitter of 6.3ms" 1 to the HIRES data. 

Again, we used our new HARPS spectra to de- 
termine wsiiij_^^2.4 ±_0.5 km s , and the calibra- 
tion of iBruntt et al.l (|2010l ) to adopt v mac — 1.0 ± 
0.3 km s -1 . This estimate of v si n I agrees well with th e 
vsinl = 2.5 ± 0.9 km s" 1 from lAnderson et al.l (|2011f ). 
but is substantially different t o the value of v s in I = 
0.4 ± 0.4kms-l obtained by IBekv et al.l (|2011h . who 
used v mac = 3.2 9 kms— 1 based on the calibration of 
iValenti fc Fischerl ((20051 ). We found that using such a 
high macroturbulence value led to a poor fit for many 
of the spectral lines, e ven with t>s i n I = 0.0 kms -1 , and 
therefore suggest that IBekv et"aTI ((201 lh have overesti- 
mated the broadening in their SOPHIE spectra. The 
IValenti fc Fischerl {2005) calibration provides only an up- 
per limit on the macroturbulence, which for cool stars 
such as WASP-40 can be significantly different from the 
true values. 

As with our RM analysis of WASP-32, we found that 
there was little to separate the different models for the 
WASP-40 system, as no significant differences were ap- 
parent in the values of Xrcd that we obtained. The ec- 
centricities returned for models with non-circular orbits 
were found to be insign ificant by the statistical test of 
ILucv fc Sweenevl (|1971[ L and the values were all found 
to be consistent with e = to within 1.5<r. We also 
note that the addition of HARPS spectrographs mea- 
surements has reduced the value of any possible eccen- 
tricity in the system by a facto r of 10 compared to the 
results in lAnderson et all (|2011| ). The addition of a long- 
term radial velocity trend to the model was found to 
provide no improvement in the quality of the fit ob- 
tained, and the low magnitude of any possible trend 
(|7| < 43ms -1 yr" 1 ) leads us to conclude that no such 
trend is present in the system. Imposing a prior on the 
stellar radius produced only small changes in the mass 
(|<5-M*| < 2 percent), radius (<5i?„ < 3 percent), and den- 
sity ((5p* < 7 percent) of the host star. The impact pa- 
rameter was similarly unaffected, with only the error bars 
increasing with the relaxation of the prior. 

Adding a prior on v sin I using the spectroscopic mea- 
surement produced no change in the value of XJ? ed , but it 
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Orbital phase Orbital phase 

Fig. 3. — : Radial velocity curve produced by our optimal model for the WASP-38 system. The model uses an eccentric 
orbit and a prior on the stellar radius, but no long-term radial velocity trend is found and the prior on the spectroscopic 
vs'ml is not applied. Data from CORALIE are denoted by triangles. Data from SOPHIE are denoted by circles. Data 
from FIES are denoted by diamonds. Data from HARPS are denoted by squares. Error bars are marked in grey; some 
are smaller than the size of the data points that they accompany. Format as for Fig.[TJ Colour versions of these figures 
are available in the online edition of the journal. 




Barycentric RV [km s 1 ] Barycentric RV [km s 1 ] 

Fig. 4. — : Left: Residual map of WASP-38 time series CCFs with the model stellar spectrum subtracted. The bright 
signature of the planet is clearly visible, and it's trajectory from bottom left to top right clearly indicates a prograde 
orbit. Right: The best- fitting model for the time- variable planet feature has been subtracted, leaving the overall 
residual map. The lack of any remaining signatures suggests that the star is chromospherically quiet. 
Details as for Fig.[2j 
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significantly lowered the value of v sin / returned by the 
MCMC algorithm, and greatly reduced the uncertainties 
on the values of A that were being produced. Exam- 
ination of the HARPS spectroscopy indicated that the 
amplitude of any RM effect was likely to be low, with 
the error bars on the data such that they obscured any 
possible anomaly in the RV curve. This indicated that 
the value of v sin / was likely to be low, and that the 
error bars on A would likely be high. This information, 
combined with the lack of any difference in the quality 
of fit, led us to select a solution in which the prior on 
vsml was not applied. 

Our adopted solution therefore uses the combina- 
tion of a circular orbit and no long-term radial ve- 
locity trend, with neither the prior on vsinl nor the 
prior on the stellar mass applied. The radial veloc- 
ity curve that results is shown in Fig.[5l Values of 
A = 24.2°+ 4 7 4 6 5 °, vsinl = 0.6±g;X kms" 1 , b = 0.87 ± 0.01 
and i — 85.0° ± 0.2 were returned for this combination 
of conditions. We note that the value we obtain for 
the impact parameter is consistent with that found by 
lAnderson et al.l (|2011| ). who found a 40.5 percent likeli- 
hood that the system is grazing. We also note that this 
system se rves as a good examp le of the systematic dis- 
cussed bv lAlbrecht et al.l (|2011l ). They showed that sys- 
tems with low-amplitude, low signal-to-noise measure- 
ments of the Rossiter-McLaughlin effect were preferen- 
tially found to be either close to aligned (A = 0°) or anti- 
aligned (A = 180°), with the posterior-probability distri- 
butions of these systems showing greater ranges of possi- 
ble solutions around these angles. Fig.[B] shows the pos- 
terior probability distribution for A against v sin / from 
our MCMC run. It is immediately clear that there are a 
greater number of solutions, covering a greater range of 
values, for v sin / at angles close to 0°; the effect at 180° is 
less pronounced. We note that our solutio n lies relatively 
close t o the former angle, as predicted bv lAlbrecht et al.l 
(2011), but we also note that our error bars are such 
that a wide range of alignment angles are included in the 
possible range of solutions that we find. 

In light of this, we analysed the system using our pre- 
ferred choice of priors and initial conditions, but with no 
Rossiter-McLaughlin effect fitting. We found that this 
produced results that showed no difference in terms of 
quality of fit from our adopted solution, with a value of 
Xred noRM = 1-3 ±0.2 that is in complete agreement with 
X?cd = 1-3 ± 0-2 from the solution adopted above. We 
therefore consider our weak constraints on the alignment 
angle to be equivalent to a non-detection of the Rossiter- 
McLaughlin effect. 



5.1. Doppler tomography 

We attempted to model the system using Doppler to- 
mography, but the combination of the low signal-to-noise 
and slow rotation proved too difficult to analyse using 
this method. This nicely highlights a major limitation of 
the technique, namely systems with poor quality spectro- 
scopic data. RM analysis is able to overcome the poor 
data quality to provide a result, although it may be in- 
conclusive. However the tomography method is simply 
unable to process the data if the effect of the planetary 
transit on the stellar line profile is insignificant. 



6. OUR RESULTS IN CONTEXT 

We now consider our new results, which are sum- 
marised in TableEl in the context of the complete set of 
spin-orbit alignment measurements. To date, 52 systems 
have such measurements published; our results push that 
number up to 54. 

6.1. Effect on existing trends with mass and 
temperature 

InlBrown et al.1 (I2012D we updated the |A| — T e // plot 
of iWinn et al.l (|2010af) with all of the systems published 
since their original analysis, and t he re sults presented 
in our own work. lAlbrecht et al.l ((2012) carried out a 
similar exercise with their new measurements. We con- 
sider the new results presented herein in the context of 
the set of systems listed in the Holt- Rossiter-McLaughlin 
database compiled by Re ne Heller 12 , as well a s the new 
systems and results from lAlbrecht et al.1 (|2012h . 

Our new results add little further information to 
the detected trend with temperature discovered by 
W inn et all (|2010af) . The three systems that we study 
herein all fit into the 'cool' category, although WASP- 
32 at T c ff = 6100 ± 100 lies close to the critical tem- 
perature dividing the two sub-populations, and WASP- 
38 at T e ff — 6110 ± 150 encompasses the value of 
T cr it — 6250 K within it's la T e g range. Similarly, our 
new results have no effect on the known trend with plan- 
etary mass as the systems all have M p < 3Mj up . They 
therefore cannot provide counter-examples, as planets in 
this category are already thought to exhibit randomly 
distributed values of A. 

6.2. Stellar ages 

The host stars of WASP-32 and WASP-40 arc in- 
sufficicntl y massi v e to f ulfil the selection criterion im- 
posed by iTriaudl (|2011l ) for his study of the trend of 
A with stellar age. WASP-38 lies close to the cut- 
off mass; in some of our simulations it fa lls below the 
limit, but in our adopted solution it fulfils iTriaudl 's cri- 
terion for inclusion. We computed the ages for our 
three systems using a simple isochrone interpolation rou- 
tine, and several different sets of stellar models, in an 
attempt to better characterise the inherent uncertain- 
ties . Specifically, we made use of the Padova mod- 
els (|Marigo et aill2008l; iGirardi et aUpOlOl) . Yonsei-Yale 
(Y 2 ) models (|Demaraue et al.l 120041 ). Teramo models 
(jPietrinferni et al.l 120041 ), and Victoria- Regina S tellar 
Structure (VRSS) models (|VandenBerg et alj|200l . We 

— 1/3 

carried out our isochrone fits in - T c s space, tak- 

ing the effective temperature from spectroscopic analy- 
sis of the HARPS spectra and the stellar density value 
as found by our preferred model under the tomographic 
method for WASP-32 and -38, and the RM model for 
WASP-40. The ages that we obtained for the three sys- 
tems using these models are set out in Table. |4j We also 
assessed the ages of the systems using a combination of 
the R' HK activity metric, and gyrochronology. We mea- 
sured the chromospheric Ca II H & K emission from the 
new HARPS spectra that we obtained for the three sys- 
tems discussed herein, using this to calculate \og(R' HK ). 
We then computed the stellar rotation period using the 

12 As of 2012 April 25. |http://www.aip.de/People/RHelTer| 
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Fig. 5. — : Results from the fit to the radial velocity data for our adopted solution for WASP-40. A circular orbit was 
used, with no prior on the spectroscopic vs'ml, no long-term radial velocity trend, no prior on the stellar radius. Data 
from HIRES are denoted by triangles. Data from SOPHIE are denoted by circles. Data from HARPS are denoted 
by squares. Error bars are marked in grey; some are smaller than the size of the data points that they accompany. 
Format as for Fig. [I] Colour versions of these figures are available in the online edition of the journal. 



method of iWatson et al.1 (|2010f ). which in turn allowed 
us to estimate the age of the system using the method of 
Barnes! (120071). coup led wit h the improved co efficients of 
Meibom et all (|2009D and lJames et all (|2010f ). 



The age estimates for the WASP-32 system vary quite 
considerably but do all agree within the la errors, with 
the age from gyrochronology lying in the middle of the 
range. However this is more a function of the rather 
large uncertainties than any indication that the age is 
well constrained. The system nicely highlights the de- 
pendence of isochronal age estimates on the set o f stella r 
models that is used. As discussed bv lSouthworthl (|2009f ). 
the choice of stellar models can have a significant impact 
on the derived properties of exoplanetary systems, par- 
ticularly through the introduction of systematic errors. 
iSouthwort h (2009) also notes though that different sets 
of models are often based on the same physical under- 
pinnings, differing only in their implementation, and that 
therefore the results canno t be considered to be totally 
independent. iTriaudl (|2011| ) pointed out that isochronal 
analysis is less precise for stars with M* < 1.2 M Q owing 
to the increased length of their main sequence lifetime, 
and their less pronounced radius increase (and therefore 
density decrease) than more massive stars. We found a 
mass of 1.11 ± 0.02 M from our RM analysis, in agree- 



ment with the value from iMaxted et all (|2010f > ; the wide 
age range that we obtain is therefore expected given our 
preceding comments. The WASP-38 system on the other 
hand shows much better agreement between the age esti- 
mates obtained using the different sets of stellar models. 
With a mass of 1.17 ± 0.02 M Q from our RM a nalysis, 
the sy stem lies closer to the arbitrary cut-off of ITriaudl 
(2011), so we might expect that the age would be bet- 
ter constrained. Nevertheless, our four age estimates for 
WASP-38 all agree with the postulated trend for align- 
ment angle to decrease with time 

Finally, WASP-40 is poorly constrained, and we are 
unable to place upper limits on the age using the avail- 
able isochrones for 3 out of the 4 model sets that we tried. 
It is hard to conclude anything from this, but the differ- 
ent models do agree that the system is older than either 
WASP-32 or WASP-38. On the other hand, we note that 
the gyrochronological estimate of the stellar age is signif- 
icantly lower than the age limit that we found from our 
isochronal fits to the Teramo and VRSS models, although 
it is con sistent with the results from the Y 2 and Padova 
models. lAnderson et all (|2011| ) found age s for the sys- 
tem o f 6 ± 5 Gyr using th e stella r models of lMarigo et all 
(|2008f ) and lBertelli et all (|2008l ). which is consistent with 
our values. They too found a lower age (l.Zt [ s Gyr) us- 
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Parameter Units 



WASP-32 3 



WASP-38 a 



WASP-40 b 



Fitted Parameters 
D 

K m s _1 

b R, 

W days 

P days 

To BJD\jtc 



0.0108 ±0.0001 
0.478 ±0.011 
0.66 ±0.02 
0.0990 ± 0.0007 
2.718661 ± 0.000002 



0.0069 ± 0.0001 
0.252 ±0.004 



0.12 



+0.08 



-0.07 

0.1969 ±0.0010 
6.87188 ± 0.00001 3.039577 



0.0143 ±0.0005 

0.0912 ± 0.002 

0.87 ±0.01 

070+ 001 
u.uiu_ Q 002 



+0.000005 
-0.000006 



2450000 5681.1945 ±0.0002 5322.1774 ± 0.0006 5407.9088 ± 0.0002 



Derived parameters 








Rp/ R* 




0.104 ±0.005 


0.083 ± 0.002 


120+ 009 
u - lzu -0 .007 


R,/a 




0.129 ±0.003 


n ns9Q+ 000S 
0.0S29_ Q 0007 


102+ 003 
u.iuz_ 004 


R* 




1.09 ±0.03 


1.35 ±0.02 


0.87 ±0.04 


M* 


Mq 


1.07 ±0.05 


1.23 ±0.04 


0.92 ± 0.06 


p* 


po 


0.84 ±0.05 


0.50 ±0.01 


-1 OO+0.16 

J-^-o.ia 


[Fe/H] 




-0.13 ±0.10 


-0.02 ± 0.07 


0.14 ±0.11 


v sin 7 


km s _1 


3 9+0.4 


' -°-0.2 


6 +0 - 7 


Rp 


^Jup,eq 


1.10 ±0.04 


1.09 ±0.02 




M P 


M Jup 


3 46+ 14 


2.71 ±0.07 


0.62 ± 0.03 


a 


AU 


0.0390 ± 0.0006 


0.0758 ± 0.0008 


0.0400 ±0.( 


i 


deg 


85.1 ±0.2 


89.5«-j 


85.0 ±0.2 


e 




0( adopted) 


0.028 ± 0.003 


O(adopted) 


W 


deg 




-22.2^-2 





X 


deg 




'■°-6.1 


24 2+ 76 ' 


\i\ 


msyr" 1 


0( adopted) 


0( adopted) 


O(adoptcd) 



"Results from Doppler Tomography. 
'Results from Rossiter-McLaughlin analysis. 



TABLE 4: Age estimates for the three systems. 



Stellar model fitting age Gyrochronology 
Padova Y 2 Teramo VRSS age 

System (Gyr) (Gyr) (Gyr) (Gyr) (Gyr) 



WASP-32 2.36tJ'^ 2.22^1 4.50t^ l-4llj;^ 2.42^°; 
WASP-38 3.41±g;J| 3.29^0.53 3.59^0.70 3.20l°:s9 3.4lt°; : 
WASP-40 > 1.20 6.36±3 11 > 4.96 > 5.73 3.60ti; 



ing gyrochronology, in their case based on an estimate 
of the rotation period derived from vsinl, but this does 
not match our estimate. From this information we ten- 
tative ly predict, following the trend noticed by iTriaudl 
(|2011h . that the system will prove to be aligned if the 
uncertainty on A is able to be reduced. 

6.3. Are the systems aligned? 

In iBrown et~aH (|20L1 we introduced a new test for 
misalignment that is based upon the Bayesian Informa- 
tion Criterion (BIC). The BIC is calculated for the set of 
radial velocity measurements that fall within the transit 
window, for both the best-fitting model and one that as- 
sumes an aligned orbit with A = 0°. The ratio, B, of the 
aligned orbit BIC to the free-A BIC is then calculated. 
Systems with B < 0.99 are classed as aligned, systems 



with B > 1.01 are classed as misaligned, and systems 
that f all between the s e limi ts are classed as indetermi- 
nate. lAlbrecht et~aTl (|2012[ ) note that the BIC test is 
affected by the relative numbers of RV measurements in 
transit compared to out of transit, and that it assumes 
that no correlated noise is present. We acknowledge that 
these are indeed shortcomings of our test, and that they 
might affect the boundari e s betw een the three categories 
discussed in IBrown et al.l ([2012), but we comment that 
the test is still quantitative, as opposed to the qualit ative 
nature of the previ ous tests in iTriaud et al.l ([2010D and 
IWinn et al.l (I2010al) . 

We applied this test to our new results, calculating val- 
ues of 0.92 for WASP-32, indicating alignment, 0.95 for 
WASP-38, indicating alignment, and 0.97 for WASP-40, 
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1 1 | II II | II II | II II | II II | II II | II II | II : 




-1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1- 


- 1 1 - 



-150 -100 -50 50 100 150 0.00 0.02 

M°) 

Fig. 6. — : Posterior probability distribution for vsinl 
and A, derived from the Markov chain, for the fit to the 
data for WASP-40 described in Fig. El The white contour 
marks the 62.87 percent confidence regions, the black, 
dashed contour the 95.45 percent confidence regions, and 
the black, dotted contour the 99.73 percent confidence 
regions. Marginalised, ID distributions are displayed in 
the side panels. A = lies well within the main body of 
the distribution. 



indicating alignment. However in I Brown et al.l (|2012f l 
we also postulated a fourth category, that of 'no detec- 
tion', defining this as vsinl consistent with to within 
lcr. With further reflection, we consider this definition 
to be inadequate. Our analysis routines are set-up in 
such a way that such a scenario is highly unlikely to 
exist; indeed, a lower error bar on vsinl of greater mag- 
nitude than the value itself is nonsensical, as negative ro- 
tation is a physical impossibility when considering only 
the magnitude of the rotation. We therefore revise this 
category of 'no detection' to include systems with vsinl 
consistent with to within 2a. This new definition en- 
compasses WASP-40, as we feel is appropriate given the 
poor signal-to-noise of the data that we obtained, and 
the indistinct Rossiter-McLaughlin effect that we find. 
The amended category include s no additional syst ems 
from our previous sample (see iBrown et al.l (|2012ft for 
details), although WASP-1 and WAS P-16 are close, but 
would include t he results for TrES -2 (jWinn et al J 12008ft 
and HAT-P-11 ([Winn et al.H2010bft . 

6.4. Tidal timescales 

lAlbrecht et alj (|2012f ) present two different approaches 
for estimating the tidal evolution timescales for hot 
Jupiter systems, and calculate said timescale for a large 
sample of planets for which the Rossiter-McLaughlin ef- 
fect has been measured. They take two approaches. In 
the first, they consider a bimodal sample of planets: 
those with convective envelopes, and those with radia- 
tive envelopes. In the second approach they consider 



the mass of the convective envelope, which they link to 
stellar effective temperature. Unfortunately this second 
approach relies on an unspecified proportionality con- 
stant, and the relation between T e g and Mcz that they 
derived is also unknown. We therefore consider their first 
approach, which is encapsulated in the equations 



1 

TCE 



1 



2 /a/i?, 



and 



1 



1 



0.25 • 5 • 10 9 yr 



10 • 10 9 yr 



40 



5/6 



a/R, 
6 



-17/2 



(2) 



(3) 



The stellar effective temperatures of our three sys- 
tems are, as mentioned previously, below the critical 
temperature dividing the 'hot' and 'cool' regions of pa- 
rameter space. They therefore all fall under the con- 
vective envelope version of the tidal timescale equation. 
Using parameters from our best-fitting models (tomo- 
graphic for WASP-32 and WASP-38, and RM for WSAP- 
40), we calculate the tidal timescales for our three sys- 
tems using H We find r CE = 5.46867636 x 10 10 yr for 
WASP-32, T CE = 1.79351519 x 10 12 yr for WASP-38, and 
T CE = 5.41028372 x 10 12 yr fo r WASP-40. These va lues 
fit nicely into the scheme that lAlbrecht et al.l ([20121 ) de- 
veloped, whereby systems in which the tidal timescale is 
short preferentially show low values of A, whereas those 
with longer timescales appear to present an almost ran- 
dom distribution of A. The timescales for our three sys- 
tems are relatively short, particularly where WASP-32 
is concerned, and the small alignment angle that we ob- 
tained for that system is exactly as expected. 

7. WHY USE DOPPLER TOMOGRAPHY? 

As we discussed in Section[TJ Doppler tomography is 
one of a number of methods for characterising spin-orbit 
alignment that are beginning to be used as alternatives 
to the traditional radial velocity based approach that 
we have used to analyse all three of the systems in this 
study. Although tomography has weaknesses, and can- 
not be applied to every planetary system (as witnessed 
with WASP-40 previously), it has one great selling point 
over the radial velocity method. Tomography is able to 
lift the strong degeneracy that exists between v sin / and 
A, and which is strongest in systems with low impact 
parameter. 

The geometry of the path that the planetary orbit 
traces across the stellar disc affects the uncertainty in 
the spin-orbit alignment angle, particularly if that path 
is symmetric with respect to the approaching and reced- 
ing hemispheres of the star. As the impact parameter in- 
creases, the range of alignment an gles than can produc e 
a symmetric RM curve decreases ()Albrecht et al.l [20 lift . 
The limiting cases illustrate this well. With 6 = 0, any 
value of A will produce equal transit path lengths through 
the red- and blue-shifted halves of the stellar disc, whilst 
with 6=1, only A = 0° and A = 180° will have the 
same effect. Thus as 6 decreases, the uncertainty in the 
estimate of A increases. 

This is not the only parameter involved however. The 
stellar rotation, vsinl, dictates the amplitude of the 
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Rossiter-McLaughlin anomaly, but this is often ambigu- 
ous owing to the uncertainties present in the radial veloc- 
ity measurements. It is often not clear, particularly for 
systems with low v sin /, whether the anomaly is truly 
asymmetric, or whether it is an effect pro duced by the 
error bars (see for example WASP-25 in iBrown et al.1 
((20121) 1. This means that the same anomaly can often 
be fit in two different ways. Either vsinl is low and A 
indicates misalignment, with the resulting asymmetry in 
the model used to fit the uncertainties, or A is low and a 
rapid v sin / is used, with the greater amplitude providing 
the required fit. Often what results is a compromise solu- 
tion, with large error bars on both parameters and some 
degree of degeneracy between them. This arises owing to 
the use of the Rossiter-McLaughlin effect to characterise 
both parameters simultaneously. The problem is exac- 
erbated for systems with low signal-to-noise, for which 
the range of possible models that fit the data is greatly 
increased owing to the greater relative size of the uncer- 
tainties, and for systems with low impact parameter, for 
the reasons discussed above. 

The Doppler tomography method does not suffer from 
this same problem, and is therefore able to provide bet- 
ter constraints on A in these problematic cases. Directly 
modelling the separate components of the CCF provides 
several separate constraints on the parameters involved 
in the model, and the geometric calculation of the posi- 
tion of the planet's shadow on the stellar disc helps to 
remove ambiguity regarding A. These two factors lift the 
degeneracy experienced with the traditional method. 

WASP-38, as an example of a system with low im- 
pact parameter, provides a reasonable example of the 
advantages that the tomographic analysis method holds 
over the standard radial velocity method. Table[2] clearly 
shows that the error bars on A have been decreased by the 
use of Doppler tomography, and Fig. [7] shows the change 
in the relationship between the values of v sin I and A for 
the two analysis methods. The two posterior probabil- 
ity density plots show completely different distributions, 
with that for the radial velocity method showing a clear 
correlation between the two parameters, with obvious 
degeneracies in the fitted values. The tomographic dis- 
tribution, on the other hand, shows very little in the way 
of correlation, and although there is still some spread in 
the A distribution the range of v sin / values has quite 
clearly been heavily restricted. 

8. CONCLUSIONS 

We have presented measurements of the sky-projected 
spin-orbit alignment angle for the hot Jupiters WASP- 
32, WASP-38, and HAT-P-27WASP-40, using both the 
Rossiter-McLaughlin effect and Doppler tomography. 
We find that WASP-32 exhibits an alignment angle of 
A = 10.5 o ^5 g (from Doppler tomography) and a ro- 
tation speed of vsinl = 3.9^q 5 kms -1 , indicating an 
aligned system. The results from our two analysis meth- 
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ods are consistent, and show good agreement, where ap- 
plicable, with the original discovery paper. For HAT-P- 
27/ WASP-40 we find a much lower rotation speed than 
suggested by the discovery paper and spectroscopic anal- 
ysis, vsinl = 0.6^q 4 kms -1 , but our poor signal-to- 
noise data allows us to place only weak constraints on 
the alignment angle. We find A = 24.2° ^l®'® , which we 
classify as a non-detection, and are unable to apply the 
tomography method to the syste m. For WASP-38 weim- 
prove on the previous analysis of iSimpson et al.l (|2011l) . 
reducing the uncertainty in A by an order of magnitude, 
and obtaining A = 7.5°^"g 1 7 and vsinl = 7.5^Q'.!,knis~ 1 
through tomographic analysis. Our results again agree 
well between the two analysis methods. 

We consider the effect of our new results on the com- 
plete sample of hot Jupiters for which A has been mea- 
sured, finding that they either provide support for, or 
no evidence in opposition to, previously existing trends 
within the ensemble. We also consider the benefits of us- 
ing Doppler tomography over the Rossiter-McLaughlin 
analysis method, and comment that its use is helpful in 
lifting degeneracies in the fitted solution that arise when 
using the more traditional method. 

NOTA BENE 

We have used the UTC time standard and Barycentric 
Julian Dates in our analysis. Our results are based on 
the equatorial solar and jovian radii, and masses taken 
from Aliens Astrophysical Quantities. 
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Fig. 7. — : Posterior probability distributions for vsinl and A for both of the analysis methods discussed in this work. 
These distributions are for the analysis of the WASP-38 system discussed in Section^ Key as for Fig. [51 Left: Radial 
velocity measurement based Rossiter-McLaughlin analysis. Right: Doppler Tomography. The difference between the 
two methods is stark, with the tomographic analysis yielding a much reduced correlation between the parameters. 
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TABLE 5: RV data for WASP-32 obtained using HARPS. 
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TABLE 6: Out-of-transit RV data for WASP-38 obtained using HARPS. 
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BJD UTC (-2450000) (kms" 1 ) (kms^ 1 ) 
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TABLE 7: In-transit RV data for WASP-38 obtained using HARPS. 
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B.TD UTC (-2450000) (kmr 1 ) (kms -1 ) 
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00618 


5727 


574140 


-9 


70082 





00615 


5727. 


578017 


-9 


.69531 





.00609 


5727 


581882 


-9 


.70552 





.00618 


5727 


585586 


-9 


.70640 





.00618 


5727 


592484 


-9 


.71415 





.00667 


5727 


596430 


-9 


.72199 





.00718 


5727 


600203 


-9 


71776 





00651 


5727 


604046 


-9 


73514 





00681 


5727 


607923 


-9 


73209 





00676 


5727. 


611800 


-9 


73769 





00745 


5727. 


619335 


-9 


74676 





00739 


5727 


623177 


-9 


75379 





.00724 


5727 


627020 


-9 


76161 





.00671 


5727 


630746 


-9 


.76295 





.00678 


5727. 


634589 


-9 


.76762 





.00661 


5727. 


638466 


-9 


.76457 





00690 


5727. 


642667 


-9 


78228 





00754 


5727 


646510 


— 9 


77732 





00655 


5727 


650422 


— 9 


78573 





00707 


5727. 


654183 


— 9 


78749 





00695 


5727. 


658164 


— 9 


78420 





00719 


5727 


661868 


-9 


79843 





00745 


5727. 


666231 


-9 


.80281 





.00769 


5727. 


670143 


-9 


.79479 





.00621 


5727. 


673905 


-9 


79231 





.00634 


5727. 


677747 


-9 


79598 





00576 


5727. 


681543 


-9 


79677 





00532 


5727. 


685386 


-9 


.80741 





00569 


5727 


689656 


-9 


.80222 





00564 


5727. 


693430 


-9 


.79689 





00612 


5727. 


697226 


-9 


79349 





00704 


5727 


701242 


-9 


.78573 





.00693 


5727 


705049 


-9 


78707 





.00585 


5727 


708915 


-9 


.78243 





.00568 


5727 


713197 


-9 


.76337 





.00552 


5727. 


716936 


— 9 


.77080 





.00593 


5727 


720871 


— 9 


77108 





00529 


5727 


724713 


-9 


76885 





00528 


5727 


728544 


-9 


76841 





00528 


5727. 


732306 


-9 


.76986 





00514 


5727 


736484 


-9 


76921 





00541 


5727 


740141 


-9 


.76496 





.00641 


5727 


744192 


-9 


.76324 





00720 


5727 


747999 


-9 


77120 





.00705 


5727. 


751981 


-9 


76748 





00679 


5727. 


755823 


-9 


76974 





.00648 


5727. 


760256 


-9 


78057 





00585 


5727. 


764376 


-9 


.78068 





00588 



The alignment of WASP-32, WASP-38 and WASP-40 



TABLE 8: RV data for WASP-40 obtained using HARPS. 



RV u RV 
BJD UTC (-2450000) (kms- 1 ) (kins- 1 ) 



5686 


.691040 


-15 


.68740 





.00531 


5691 


.594153 


-15 


.83574 





.00593 


5691 


621178 


-15 


.84591 





.00639 


5691. 


.784798 


-15 


.82559 





.00548 


5692. 


.677700 


-15 


.67883 





.00455 


5692. 


.747097 


-15 


.67360 





.00411 


5693 


.571851 


-15 


.75439 





.00814 


5693 


.577950 


-15 


.75479 





.00878 


5693 


.584107 


-15 


.74939 





.00838 


5693 


.590207 


-15 


.75567 





.00843 


5693. 


.596480 


-15 


.74962 





.00807 


5693. 


.602579 


-15 


.75667 





.00803 


5693. 


.608737 


-15 


.76548 





.00767 


5693 


.614952 


-15 


.75844 





.00817 


5693 


621167 


-15 


.76064 





.00735 


5693 


.627093 


-15 


.75956 





.00822 


5693 


.633296 


-15 


.77570 





.00901 


5693. 


.639569 


-15 


.76272 





.00858 


5693. 


.645669 


-15 


.77900 





.00757 


5693. 


.651884 


-15 


.77096 





.00700 


5693. 


.657925 


-15 


.77981 





.00764 


5693 


.664152 


-15 


.76408 





.00737 


5693. 


.670263 


-15 


.76868 





.00751 


5693 


.676467 


-15 


.76874 





.00774 


5693. 


.682682 


-15 


.76226 





.00709 


5693. 


.688781 


-15 


.78715 





.00662 


ouyo. 
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5693. 


.701142 


-15 


.78574 





.00677 


5693 


.707184 


-15 


.78466 





.00663 


5693. 


.713387 


-15 


.76800 





.00666 


5693 


.719603 


-15 


.78821 





.00713 


5693. 


.725633 


-15 


.77244 





.00792 


5693. 


.787541 


-15 


.78723 





.00474 


5694. 


.567017 


-15 


.85210 





.00448 


5695 


.576799 


-15 


.69093 





.00362 


5695. 


.757478 


-15 


.67832 





.00406 



